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I PROGRESS OF RESEARCH 

A. Introduction 

This semi-annual report describes continued progress on the research. Among several 
approaches in this area of research, our approach to the parametric uncertainties 
are being matured everyday. This approach deals with real parameter uncertainties 
which other techniques such as H°° optimal control, fi analysis and synthesis, and 
t} optimal control cannot deal. The primary assumption of this approach is that 
the mathematical models are well obtained so that the most of system uncertainties 
can be translated into parameter uncertainties of their linear system representations. 
These uncertainties may be due to modeling, nonlinearity of the physical system, 
some time-varying parameters, etc. 

In this report period of research, we are concentrating on implementing a com- 
puter aided analysis and design tool based on new results on parametric robust sta- 
bility. This implementation will help us to reveal further details in this approach. 

B. Computer Aided Analysis and Design: Parametric Robust Stability 

There are two basic available frameworks: coefficients of the transfer function and 
parameters of the transfer function. If coefficients of the characteristic polynomial 
contain parameters of interest and these coefficients are subject to perturbations 
independently, we use the coefficient perturbation framework. However, this situation 
is very rare and not very realistic. More general setup is the case of coefficients being 
linear combinations of parameters of interests. This formulation fits the case of SISO, 
SIMO and MISO. Of course, in general a MIMO system provides the characteristic 
polynomial whose coefficients are nonlinear functions of parameters. So far there is 
no available result to directly handle this case. However, by accepting a reasonable 
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amount of conservatism, one can easily reduce this situation to the second case we 
described above. 

In this section some examples are given which demonstrate how this computer 
aided tool is working. The package is incorporated with the well known software 
MATLAB in order to improve portability. The package has not yet been completed, 
and is still under development. 


1. Interlacing Property of Single Polynomial 
Consider the following polynomial, 


6(s) : = So + 8 x s + S 2 s 2 -f • • • + S n -\8 n 1 + S n s n 

= So 4" S 2 s 2 4" "4- • • • 4* S x s 4" ^s 3 4" S^s^ -\- 

' v ' ' v 

Severn (^) 

: = 6(s) \ l=ju 


6odd (*) 


— S 0 — S 2 u 2 + S 4 oj 4 — Yjuj (S x — S 2 iv 2 4- — •••). 

' V ' ' V ' 

Theorem 1 The polynomial 6(s) is Hurwitz stable if and only if 


a) ^cven(^) arccJ ^oda(^) have only simple roots and these roots interlace, 
h) For all u € 71, ^odd( a, )^«v«( w ) - ^odd^nC^) > 

The figure 1 shows the interlacing property of Hurwitz polynomial 

8(s) = 6 + 49a 4- 155s 2 4- 280s 3 + 331s 4 + 266s B 4- 145s 6 4- 52s 7 + 11s 8 + s 9 


Next some of the values of the coefficients are increased to observe how the graph 
changes. We selected and ^8- The figure 2 shows the changes in the graph 

and the given polynomial becomes unstable when it violates the interlacing property. 
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2. Interlacing Property of Interval Polynomial 

The well known Kharitonov’s theorem[l] can be easily verified by the interlacing 
property. Consider the following 

5(a) := 5 0 + S \3 -f 5 2 a 2 + • • • + 5 n _js n 1 + S n 3 n 


where 


Si < Si < Si 


and its four Kharitonov’s polynomials 

AT^a) := So + Si 3 + S 2 s 2 + 5 3 a 3 + ^a 4 + fca 5 + 5 6 a 6 + • • • 

= Sq + S 2 s 2 + ^a 4 + Ses 6 + • ■ • + S 1 3 + 5 3 a 3 4- S 5 33 4 S?s 7 -f 

' v ' ' v 

K 2 (s ) So 4 Sis + S 2 s 2 + S 3 33 + Sa s * 4" 5 B a 5 + Ses 6 + • ■ • 


K 3 (s) := 


K 4 (s) := 


So 4 " S 2 s 2 + £4 + 5 B a 6 4 • ■ ■ 4* 5ia + 53 a 3 4 5 B a 5 + S 7 37 4 • ■ 

*asiw ww 

5o 4 S 1 3 4 S 2 3<i 4 S 3 s 2 4 54 a 4 4 ^a 6 4 Se s ° 4 • • • 

S 0 4 5 2 a 2 4 64 a 4 4 S 3 36 4 • • • 4 Si 8 4 5 3 a 3 4 5ga 5 + 67 a 7 4 

' v ' ' V 

S 0 4 S\3 4 5 2 a 2 4* S 3 33 4 S 4 s* 4* 5 B a 5 4 Se 3 ^ 4" • • • 


= Sq 4" 5 2 a 2 4 S 4 s* 4 5ga 8 4 • • • + S\s 4 S 3 33 4 5 B a 5 4 Sjs 7 4 




*ssrw 


The figure 3 shows each tube constructed by K£££(u>) and K™£(u>) and 

respectively. According to the interlacing property, in order to ensure the 
stability of the family of the given polynomials, two tubes must interlace. 

The figure 4 shows the development of instability in some members in the family 
when selected coefficients are increased. 


These can also be shown by plotting perturbation boxes in the complex plane. 
The figure 5 show the traces of boxes while u; moves 0 to 00 . Each box has its vertices 
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at K%£(u>), l (u>), Kl dd(oj) and AT££d(u>). If any box contains the origin, it means 

that there is at least one unstable polynomial in the family. The figure 5 shows the 
Hurwitz stable polynomial and the figure 6 shows that the family becomes unstable 
when certain coefficients are increased. 

3. l 2 Stability Margin in Coefficient Space 

For the given Hurwitz polynomial, one often wants to know how much coefficient 
perturbation can be allowed while the family of polynomials maintains Hurwitz 
stability [2], This can conveniently be measured in terms of l 2 norm such as: 

P'= ||A*||, 

where 

r 1 r 

AS := A8\ AS 2 ■ • ■ • ■ • A8 n 

and 

8(s, AS;) := (So ± ASo) + (£1 ± Atfi)s + • • • + (£« i A<f>„)s n 

For the case of nonmonic polynomials (i.e., AS„ ^ 0), the l 2 stability margin is 
given 

p = min{6 0 > S n , 8„}. 

The expression of 8 U is found in [2]. The figure 7 shows the graph of 8 U and * indicates 
the minimum value for our example 

8(s) = 6 + 49s + 155s 2 + 280s 3 + 331s 4 + 266s B + 145s 6 + 52s 7 + 11s 8 + s e 
If we consider the monic polynomial (i.e., 8 n = 1, A£ n = 0), we have 

p = min{£ 0) ^tu} 

Again, the expression of 8 U for the monic case is found in [2]. 
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4. Generalized Kharitonov’s Theorem 

Kharitonov’s theorem is powerful and elegant, but it is not very useful for studying 
control systems because it assumes that all coefficients of the characteristic polynomial 
perturb independently. An improved version of this theorem was introduced in [3]. 
For the illustrative purpose, we give the simplest version. 

For the given interval polynomials 

Pi(s) = al + aj|s -|- a\s 2 + ajs 3 + fa^j" 

P 2 (s) = a„ + ajs + als 2 + a^s 3 H (- 

where 

&i < a* < a?, i = = 1,2 

and the give fixed polynomials 

Qi(s) = rt + /3}s + /3y+fts 3 + ---+j9 1 qi s« 

Q »(•) = A a + /J?- + AV + /9g-» + ... + *-• 

The problem is to check the stability of the family of polynomials 

Qi(s)Pi(s) + Q 2 (-s)F 2 (s). 

Let us first define 

Ki(>) = *Ln,™(’) + 

«’,(•) = 

Kiw = k 
k;w = 


Step 1: Set i\(s) = Kl(s). 



Step 2: check the stability of each following segment. 


s; 1 = Q,(>)P,(s) + - A + AJTft,)] 

S = <J,Wi>,(*) + <J,(«)I(1 - A)AT?(») + Aff|(,)] 

SU = QMPM + «.(-)[(! - A)Xj(.) + AJT*(«)] 
s;, = QMP^S) + <?,(.)[( 1 - A)JT|(») + AtfJ(j)] 

s, 1 , = «i(*)Pi(») + 9»WI(1 - A)jr;w + a*J(.)] 

= Q,(s)P,(>) + <? 2 «A?W - \Qi(s)K?(,) + \Q 2 (s)Kl( a ) 

= (1 - A) [«,(.)?,(») + Q,{s)Kl(s) I +A [<3,(«)P,(a) + Q 2 (s)K\(,)} 

V . , H ■ -■ V ^ — ' 

£i(*) $a(*) 

then call the “segment lemma” with <5i(s) and ^ 2 (^)- 

Step 3: Set P-i(s) = repeat Step 2 for Sl k with k = 1,2, 3, 4. 

Step 4: Set P\(s) = Kl(s), repeat Step 2 for S with k = 1,2, 3,4. 

Step 5: Set Pi(^) = /^(.s), repeat Step 2 for S with k = 1,2, 3, 4. 

Step 6: Set P 2 (s) = K\ ( 5 ). 

Step 7: check the stabihty of each following segment. 

S?i = Gi«[(l - A)JTi 1 (j) + AJiTjW] + Q 2 (s)P 2 {s) 

S 2 n = Q 1 (s)[(l-X)K 1 1 (s) + XKl(3)] + Q 2 (s)P 2 (s) 

Sh = Qi(s)[(l-X)Kl(s) + XKl(s)] + Q 2 (s)P 2 (s) 

S 2 u = Qi(sm-X)Kl(s) + XKl( S )] + Q 2 (a)P 2 (s) 

Step 8: Set P 2 (s) = K%(s), repeat Step 7 for S% k with k — 1,2, 3,4. 



7 


Step 9: Set P 2 (s) = K 2 (s), repeat Step 7 for with k = 1,2, 3, 4. 

Step 10: Set P 2 (s ) = Kl(s), repeat Step 7 for S\ k with k = 1,2, 3, 4. 

This algorithm can be easily extended to the case of X^«Qi(' s )^ > i( 3 )- The following 
figures show the family of polynomials as functions of c o. The entire family is Hurwitz 
stable if and only if the family does not contain the origin for all ui. Clearly Figure 
9 shows the family is stable and Figure 10 shows some members in the family are 
unstable. 


5. l 2 Stability Margin in Parameter Space 

Consider the following polynomial with parameters p 1} p 2 , ■ • ■, pi 


P) : — [«oi(Pi + Api) + a 0 2(p2 + Ap 2 ) + • ’ ■ + &oi(Pi + Apj)] "b 

[ a n(Pi + A p^ + a 12 (p 2 + Ap 2 ) + h ou(p, + Ap,)]s + 

■ • • + [a„i(pi + Api) + a n2 (p 2 + Ap 2 ) + h a^fp, + Ap,)]s n 

= ( a 0lPl + O02P2 + 1- O.01PI ) + • ■ * fanlPl + <*n2P2 + h an/Pi) 3 ” 

' V ' 

s°(*) 

l 

+ Xy ( a ° fe "b + • ” + °n*-S n ) A Pk 

*=l' 77 ' 

»*(») 

' -v ' 

Af(#,Ap) 

Then the following algorithm provides the i 2 stability margin in the parameter space. 


a k (juj) = a ok + jai k u - a 2k u> 2 - ja^u 3 + a ik u 4 + ja 6k u ; 5 - a 6fc o; 6 - 


Let 




a lT (ju) a 2r (ju) a 2r (juj) 
au{j u ) a 2i{jv) a 3 ,(ju;) 


a lr (ju) 

an(juj) 


3 , 5 

1 + a 5k u; ~ 

- ajkw 7 H ) 

-V- 

aki(jv) 


b(ju>) := 

-6°{ju) 

-6°U U ) 


A(ju) := 
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Ji@i + Pi($i — ^ 2 ) + P2(@l — $ 2 ) — T c 
J282 + Pl(^2 — $l) + P2^2 — $l) = T c 


G(s) = 

s 2 {s 2 + 2p x s + 2 p 2 ) 
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0.18 <p 2 < 0.3 

0 . 04 ^ < Pl < 0 . 2^1 

In the following figure, the dotted box indicates the range of parameter perturba- 
tions to be tolerated. The circles indicate the l 2 stability margin of the characteristic 
polynomial. If a circle completely covers the box, the stability of the closed loop 
system is guaranteed under the given parameter perturbations. Figures show that 
several circles that correspond to different controllers. 



C. Appendix 
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Polynomial is Unstable 



Figure 2 
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II DISCUSSION AND DIRECTION OF RESEARCH 
The computer aided tool currently under implementation is very useful for future 
research on this topic. In addition, it is expected that the tool will be valuable for 
engineers who actually perform design and analysis of systems. 

Different aspects of robust control research are also under study. While dealing 
with various forms of system uncertainties is important, the problems currently under 
study are also important and meaningful in practical control systems. Two problems 
currently under investigation are, zero assignment and LQG/LTR and the sensor 
failure problem . Details of these problems will be discussed later. 



